ON THE STAR FORMATION RATE AND TURBULENT 
DISSIPATION IN THE GALACTIC MODELS 
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We suggest a model for star formation function and a model for dissipation of the 
turbulent energy of interstellar medium. Star formation function takes into account 
the effect of turbulization of the ISM. It is shown that application of mentioned 
relations to the hierarchical scenario of formation of galaxies allows to explain the 
observed delay of star formation in the Galaxy that corresponds to the range of 
stellar ages from 8 — 9 to 10 — 12 Gyr. 
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I. INTRODUCTION 

Observations of the galaxies reveal a great 
variety of physical processes occurring in 
them. However, it is possible to construct 
models of the galaxies that contain fairly 
small number of parameters. For instance, 
in model [1] evolution of interstellar medium 
(ISM) is defined by the rate of supernovae 
which depends on star formation history as 
well a by dissipation of turbulent energy. 

The models of galaxies with explicit mod- 
eling of ISM by gas-dynamical methods, al- 
low to resolve the regions with the scale not 
exceeding 10 pc, while star formation occurs 
over much smaller scales. Complex structure 
of the ISM and a wide range of tempera- 
ture and density in the latter allow only phe- 
nomenological approach to the description of 
star formation in the models that use explicit 
modeling. Over the time, several functions 
that determine star formation rate (SFR) 
were suggested (see review [2]). As an ex- 
ample we can list several star-formation laws 
similar to the observations-based Schmidt 
law [3] 

SFR oc . (1) 

The first example is the Kennicutt law [4] 
which is often used in the numerical models 



of galactic evolution [5, 6]: 

SFR oc ^ oc p3/2 



(2) 



The second example is dependence that is 
used in the models of evolutionary synthe- 
sis [2]: 



SFR oc e-*/"=f 



(3) 



The third example is given by the model of 
star formation which is controlled by ioniza- 
tion balance [7]: 



SFR oc p2 



(4) 
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As a shortcoming of Schmidt-type laws one 
may consider the fact that in these models 
the rate of transformation of gas into stars 
depends on the local density of the gas only 
and the effect of supernovae explosions that 
enhance the turbulent energy of ISM is lost; 
this results in decrease of the SFR. 

Dissipation in the ISM is taken into ac- 
count by consideration of shock waves [8] , in- 
clusion of artificial viscosity [9, 10] and radia- 
tive cooling [11]. However, in the modeling of 
these effects the structure of the ISM is not 
taken into account. 

In Section 2 we analyze conditions for star 
formation in the galaxies, derive a model for 
star formation and dissipation of the tur- 
bulent energy of the gas, and describe a 
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single-zone model for the evolution of galax- 
ies and its generalized version. In Section 3 
we present results of application of suggested 
SFR-function to the single-zone model of the 
Galaxy. 

II. EVOLUTIONARY MODEL OF A 
GALAXY 

In this Section we construct a model for 
galactic evolution based on the Tutukov- 
Firmani single-zone model [1, 12]. In the new 
model we use a more accurate star formation 
function and apply a mechanism for dissipa- 
tion of turbulent energy that takes into ac- 
count the structure of ISM; as well, we take 
into account viscosity. 

A. Star formation rate 

Since one of the aims of the paper is 
derivation of the star formation model for ap- 
plication in the numerical models, it is rea- 
sonable to consider ISM over the scale of 
the order of the minimum scale resolved in 
the numerical models of the galaxies - 10- 
— 100 pc, and over a time span that allows 
to consider ISM as stationary. 

Let consider SFR in the form 

i> = c^— , (5) 

Tflf 

where p is gas density, rgf is the free-fall time; 
c* is so called dimensionless efficiency of star 
formation [6]. The values of this constant 
that are encountered in the galactic models 
vary from 0.1 [13] to 1 [6]. Since it is deemed 
that the rate of star formation is proportional 
to the amount of cold molecular gas [14], it 
is reasonable to define the dimensionless effi- 
ciency of star formation as a fraction of gas 
that is cold and dense. 

Let assume that the density perturbations 
in the gas are distributed in the mass in- 
terval (mmin, ''^max) and follow a power law, 
i.e. P{dm} oc dm^~^. Then the fraction of 



mass contained in all perturbations that have 
Jeans mass or higher one, may be obtained by 
integration of mass spectrum over (mj, mmax) 
range: 

^ _ 1 - (mj/mmax)^~^ / rrimin V"^ 

1 - (mmin/mmax)^"^ \ rUj J 

(6) 

If mmin and mmax ai'c of the order of minimum 
and maximum stellar mass respectively, i.e. 
0.1 — IOOM0, one may assume mmin/^^max 
0; then 

oc mj~^ . (7) 

Inserting into Eq. (5) expressions for Jeans 
mass mj oc p~^l'^T'^l'^ [15] and free-fall time 
Tff oc [6], one obtains: 

^ = ^pi/3+iri(^-^) , (8) 

where g — is a normalization constant (which 
will be defined below) and T is temperature. 

It is necessary to clarify our notion of the 
temperature. Virial internal energy of the in- 
terstellar gas per particle, averaged over all 
components of the present-day ISM is ap- 
proximately 4 X IO^/cb (where is Boltz- 
mann constant). This energy includes both 
thermal and turbulent energy. In the star for- 
mation region directly the temperature is of 
the order of tens or several Kelvins [16]. Since 
dependence (10) is suggested for inclusion in 
numerical model, we have to proceed from 
the fact that at hand there are only quanti- 
ties averaged over the length scale of compu- 
tational cell of 10—100 pc, despite star forma- 
tion occurs over much smaller length scales. 
However, assuming power-law for the mass 
spectrum of ISM components, as above, as 
well as power-law dependence of dispersion 
of turbulent velocity on the length scale [17], 
one may claim that the average the turbu- 
lent energy of the ISM is proportional to the 
turbulent energy and temperature over small 
length scales, i.e., in the star formation re- 
gions. Below, the temperature T, will imply 
a quantity that is related to the maximum 
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value of dispersion of the turbulent velocity 

2 



(75 as 



kBT 
1^ 



(9) 



where /i is the average molecular weight. 

If the power of mass spectrum is Salpeter's 
one {l3 = 2.35), ^ = ^p2.i752--2.025_ g^^Q^ 

will assume for the ISM the dependence 



(10) 



We should note that quadratic depen- 
dence of SFR on density was assumed be- 
fore in the model of star formation governed 
by ionization [7] and it recommended itself 
well in the single-zone model of galactic evo- 
lution [12]: 

V'^/p', / = 2 X 10^ cm=^ g-i s-^ (11) 

Proportionality constant g may be found 
from condition of the equality of SFR in the 
models (10) and (11) at virial temperature: 

^(4 X 10^)-^ = / , (12) 

which gives 

^ = 3.2 X 10^^ cm^ g-^ s"^ . (13) 

Using Jeans criterion in the model for 
SFR, one may generahze the star formation 

law to the case when account of the chemical 
composition of the gas (its molecular weight), 
strength of magnetic field or rotation of a 
galaxy may be important. 



This approach assumes that interstellar gas 
forms a fractal structure. Formation of 
molecular clouds in the framework of frac- 
tal approach may be modeled as aggregation 
of mini- clumps of gas outflowing from red gi- 
ants [21]. 

Within magnetohydrodynamical ap- 
proach, the structure of ISM is governed by 
magnetic field. The key role in formation 
and evolution of giant molecular clouds 
is accomplished by instabilities (thermal, 
Parker ones), while the formation of the 
nuclei of molecular clouds is controlled by 
ambipolar diffusion [22-24]. 

In the graviturbulent model, ISM is a com- 
plex of random gas flows over the length 
scale of the order of Galactic disk thickness. 
The energy is transferred from large scales to 
small ones; moreover, gas dynamics at large 
length scales is defined by supersonic turbu- 
lence and over small scales it is defined by 
gravitation [24-25]. 

Magnetohydrodynamical and graviturbu- 
lent models are the most elaborated ones. 
While all three models have their own ad- 
vantages and shortcomings, we shall use tur- 
bulent model as the base. 

Our derivation of the rate of dissipation of 
turbulent energy in the galactic gas will be 
based on the assumption that the turbulence 
exists within certain range of length scales 
(^min, ^max) and that turbulent medium may 
be represented by an ensemble of the clouds 
with mass distribution [24]: 

dM^"" 

P{dM} = — , a ^ 1.5 (14) 



B. Dissipation of turbulent energy 

Interstellar medium is a complex structure 
with energy, ionization, and gas density vary- 
ing over wide range. One may note three ap- 
proaches to the description of the processes of 
formation and evolution of ISM: "graviturbu- 
lent" , magnetohydrodynamical, and fractal. 

The fractal approach is elaborated in [18- 
20] and other studies (see references in [21]). 



and with power-law dependence of the aver- 
age density of the cloud on its size [24] : 



Pi = Po 



1.1 . 



(15) 



Let us also assume that the dependence of 
the dispersion of turbulent velocity computed 
over a volume with length scale / is a power- 
law function, as it is confirmed by numerical 
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models [17, 26] and by observations [16]: 



I 



p fti 1 . 



(16) 



In these formulas, the quantities po and ctq 
may be considered as averages over a volume 
with length scale /max- 

Let assume that the dissipation of turbu- 
lent energy occurs via collisions of clouds, 
their consequent compression by passing 
shock waves, and radiation of thermal energy. 
It is evident that the efficiency of radiation 
and, hence, the efficiency of dissipation will 
be defined by the ratio of the duration of col- 
lision and the time of radiation of thermal 
energy, i.e. cooling time [27]. If it is assumed 
that the velocity of the size / clouds is related 
to the dispersion of velocities as vf = a^ — af, 
the volume density of turbulent energy asso- 
ciated with the clouds of the size in the range 

{I, y) is 



Po- 



vf + a 



^P{d/} = Po^P{d/} 



(17) 



where P{d/} is distribution of the clouds over 
the size. The fraction of energy that is radi- 
ated during collision that lasts for Tcou,/, if the 
latter is small, is equal to the ratio of coUision 
duration and the cooling time: 



Pi'yll'i 



(18) 



where A is cooling function per pair collision, 
11 is molecular weight of the gas, psh,/ = ^iPi is 
the density of the cloud after collision (after 
the passage of the shock wave), while coef- 
ficient defines the density jump at shock 
wave. Duration of collision Tcou,; — l/Di de- 
pends on the speed of the shock wave Di 
which may be expressed via velocity of the 
cloud as Di = TjiVi. In the general case, the 
rate of radiation of turbulent energy at all 
length scales is equal to 



(19) 



where ra = A/3/(27rGpo) is the timescale of 
intervals between cloud collisions [1, 28]. 

The spectrum of clouds dimensions P{dl} 
may be found from the following considera- 
tions. Let represent the distribution of clouds 
over mass (14) as a combination of condi- 
tional probability P{dM|/} and distribution 
of the clouds over their size: 



P{dM} 



/■(ma: 

'min 



P{dM\l}P{dl} , (20) 



and let write down the distribution of the 
clouds over mass as 

P{dM|/} = 6{M - pil'^) dM . (21) 

Then, by virtue of (14) and (15) we obtain 

d/i-^ 



P{dl} = 

(a-l)(3 



Jl-X 
max 



mm 



(22) 



1 ^ 1.95, M„ 



where A 

PO^max' -^^min -^^max(^min/^max) 

Thus, the expression for efficiency of dis- 



sipation Cd that enters (19) as Q 
will be 

1-A 



/Wo 
2th' 



/1-A _ /l-^ 
''max mill 



dll 



1 — exp 



(23) 



For a mono-atomic gas the value of does 
not exceed 4, for diatomic one it does not 
exceed 6 [29]. These limits are attained in 
strong shock waves only. The velocity of 
shock wave Di may exceed the velocity of 
colliding clouds by one to one and a half 
orders of magnitude only, therefore we may 
assume that approximately 2^f/rii fa 1. It 
appeared that the spectrum of clouds sizes 
P{d/} has a strong maximum at I = /min! 
besides, /mm/^max -C 1 and with sufficient ac- 
curacy we may use for the spectrum P{d/} = 
S{1 — Zmin) dl. Finally we get 



ed = 1 — exp 



Po 



P' 



(24) 
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This function is obtained presuming that col- 
hding clouds are of the same kind, i.e., they 
have equal size. If colliding clouds have differ- 
ent size, the strong, exponential, dependence 
on arguments will relax. 

Expression (24) contains two important 
parameters: cooling function A and maxi- 
mum length scale of turbulence /max- Cool- 
ing function for the gas with the temperature 
> 10^ K weakly depends on the latter and 
it may be taken as 10^^^ erg cm^'^ s^"*^ [30]. 
The length scale of turbulence was defined 
by trial and error method and it is equal to 
5 pc. One may justify the order of magni- 
tude of the latter scale by the fact that it has 
to be limited by the size of stellar structures 
that form in the ISM, for instance, open clus- 
ters which have typical length scale of several 
pc [31]. Specific value of the maximum turbu- 
lence length scale was defined from the con- 
dition of coincidence of the time span of ob- 
served and modeled interruption of star for- 
mation process (see Sect. 3). 



C. Single-zone model 

In the Tutukov-Firmani single-zone model 
of galactic evolution [1, 12] the mass of the 
gaseous disk of the galaxy is defined by star 
formation process (which has rate ijj), return 
of gas to the ISM by evolved stars, accre- 
tion of the intergalactic gas (dM"Ydt), and 
by the mass loss via stellar wind, as well as 
by sweeping of the dust by stellar radiative 
pressure and galactic wind (dM°"*/dt): 



dM 



dt 



= --0 + 



dm0(m)[m — mi.(m)] 



X ^[t - T{m)] + j;— . (25) 



dt 



dt 



Equation for the energy balance takes into 
account energy input from supernovae and 
dissipation of energy in clouds collisions: 

—rr + ~rr ~ ^sn-E^sn-Rsn , (26) 

dt dt Td 



where K and W are total kinetic energy 
and gravitational energy of the gas, respec- 
tively; esN is the fraction of supernova en- 
ergy (-Esn) transferred to the gas per su- 
pernova (it is assumed in the model that 
esN^^SN = 0.05 X W erg). 



10"^ 

RSN = i^i't - '^SNl) 







mmax 



+ / dm(j){m)ilj[t -T{m)] (27) 

'miow 

is the rate of supernovae, tsni = 10^ yr is the 
delay of type la supernovae [32] , 



(28) 



is the timescale of dissipation of the turbulent 
energy in clouds collisions [1, 28]. 

As expression for gravitational energy, Tu- 
tukov and Firmani [1] suggested to use 



W 



GMgMH^ 



(29) 



Here Mq is the total mass of the galaxy, Hg 
is the thickness of the galactic disk. In this 
model it is assumed that the galaxy is per- 
manently in the state of virial equilibrium: 



K = -\W\. 
2' ' 



(30) 



The virial relation allows to obtain an equa- 
tion for the thickness of the disk: 



dH _ esN-E'sN-RsN-R^-f^s 
Id ~ 



5GMgMH 



3H 
10^ 



(31) 



A consequence of the condition of virial equi- 
librium is that the state of the system is de- 
fined by two quantities only - by the mass 
and thickness of the gaseous disk. 

In the model, abundances of different 
chemical species are computed too. For a de- 
tailed description of the model see [12]. The 
model has a minimum set of free parameters 
and allows to follow all integral parameters of 
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the galaxies, such as SFR, himinosity, metal- 
hcity. According to the computations carried 
out in [12, 33-35], Tutukov-Firmani model 
reproduces correctly the star formation his- 
tory, the history of chemical enrichment of 
the galaxies and of the intergalactic medium. 

A peculiarity of Tutukov-Firmani model is 
that the density of gas and its kinetic energy 
are not independent parameters, by virtue 
of virial relation (30). It is easy to modify 
the model in order to remove this peculiar- 
ity. With this aim, one may introduce an 
additional variable - turbulent energy of the 
gas U and to use Newton [36] equations to 
describe the motion of the disk boundary. In 
the Newton equation, in addition to gravita- 
tion and pressure caused by inner gas energy, 
it is necessary to take into account gas vis- 
cosity: 



dH 



V, 



(32) 



'dt ~ M[H~~dH 



V 



' (33) 



where the first term contains the force which 
is equivalent to the pressure force and grav- 
ity, V is viscosity coefficient, the last term 
corresponds to the Universe expansion in the 
standard cosmological model [37]. Conserva- 
tion law for turbulent energy appears in the 
form 

dU U MV^ ^ ^ 

-T— = —U 77 — Cd \- ^ r,TT9 + ^SN-C/SN-RSN ■ 

dt H Td 2H^ 

(34) 

In this equation the term associated with dis- 
sipation (the second term) is written down 
taking into account the efficiency of dissipa- 
tion for which an expression was derived 
in the previous Section. We should note that 
the galactic evolution model that takes into 
account dynamic terms similar to (33) and 
(34) but does not account for viscosity, was 
already derived in [38]. 

It is also possible to write down somewhat 
more accurately the term for gravitational 



force that acts upon gaseous disk from disk- 
like distribution of the mass in the galaxy: 



dW 
'dH 



2GM H 
lis £1, 



i?2 



dm . 



X 



1 - 



H 



where Mg, Hg are, respectively, the mass and 
semi-thickness of the stellar disk, Mdm, ^^dm 
are the mass and the length scale of the dark 
halo distribution. 

Viscosity coefficient u may be defined as 
a square of characteristic gas velocity in the 
Galaxy, multiplied by the free-fall time: 

1/ ^ (100 km s~^)Vff 

- 10^^ - 10^° cm^ s"^ . (35) 

We assumed u — 10^° cm^ s"^ in our compu- 
tations. 

For the star formation law (10) it is neces- 
sary to know the temperature of the gas. Let 
relate turbulent energy and mean tempera- 
ture in the gas as (see (9)) 



2 II 



(36) 



III. FORMATION OF THE STELLAR 
POPULATION OF THE GALAXIES 

Observations of our and external galax- 
ies reveal different populations of stars: halo 
stars, stars of the thick and thin disks and 
bulge. Galactic evolution model should de- 
scribe the origin and properties of differ- 
ent populations in a uniform way. Mod- 
ern models of galactic evolution fall into two 
classes [39]: the first class contains the mod- 
els of monolithic collapse with smooth evo- 
lutionary transition between stellar popula- 
tions, while into the second class belong mod- 
els in which a galaxy forms by accretion 
of separate fragments that experienced an 
epoch of independent evolution and disk for- 
mation in them is a side effect of accretion. 
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FIG. 1: Evolution of the Galaxy in the scenario of monolithic collapse (model A); t - age of the 
Universe; SFR - star formation rate; H - semi-thickness of the disk; M/Mq - the ratio of gas mass 

and total mass of the galaxy; T - temperature of the gas defined by formula (37); V - velocity of 
the disk border, ed - efficiency of dissipation. Dashed line shows results of computations for the 
standard model [12]. 



The hierarchical model of formation of galax- 
ies belongs to the latter class of models. 

The history of star formation in the galax- 
ies may be strongly non-monotonous and to 
reveal several star formation bursts in the 
lifetime of a galaxy. Star formation bursts 
may be caused, for instance, by gas accre- 
tion from intergalactic medium, by mergers 
or close passages of the galaxies. Supernovae 
outbursts may delay star formation [40]. In 
the star formation history of our Galaxy one 
may discover a pause which is seen from 
the distribution [Fe/O] [39] and from depen- 



dencies [Mg/Fe] - [Fe/H] [41] and [Eu/Ba]- 
— [Fe/H] [42]. This pause corresponds to the 
range of stellar ages of 8 — 9 to 10 — 12 Gyr 
and may be interpreted as an interruption of 
the star formation process between the end 
of formation of thick disk and the beginning 
of formation of thin disk [41, 42]. 

The mechanism for cessation of star for- 
mation may be exemplified using a closed 
(without matter exchange with inter-galactic 
medium) single-zone model, appended by 
equations (32) and (33), henceforth - model 
A. Figure 1 shows results of computations for 



8 




2 2.5 3 3.5 2 2.5 3 3.5 -1 -0.5 0.5 1 

t, Gyr t, Gyr [Fe/O] 

FIG. 2: Evolution of a satellite galaxy in the scenario of monolithic collapse (model B); Z — 
metallicity; the rest of notation as in Fig. 1. Dashed line shows results of computations for the 
standard model [12]. 
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FIG. 3: Same as Fig. 2 for the evolution of the Galaxy in the hierarchical scenario (model C). 



the scenario of monolithic collapse with the 
following initial conditions: galactic radius - 
20 pc, semi-thickness of protogalactic cloud - 



2 X 10^^ Mq, mass of the dark matter over the 
length scale of galactic disk is zero, gas tem- 
perature - 4 X 10^ K, the age of the Universe 
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that corresponds to the beginning of compu- 
tations - 3 Gyr, {z = 1.8 in the standard 
cosmology [37]), the rest of the parameters 
correspond to standard model [12]. After the 
phase of collapse which results in the increase 
of gas temperature by an order of magni- 
tude and a half, star formation interrupts for 
about 1 Gyr and then a burst of star forma- 
tion occurs with a consequent decay. As the 
plots [0/H]-t, [Fe/H]-t, and [0/Fe]-[Fe/H] 
show, evolution of chemical abundances also 
reflects the presence of the delay of star for- 
mation. However, observations suggest that 
iron abundance that corresponds to star for- 
mation delay is [Fe/H] > —0.5 [41], while in 
the model under consideration cessation of 
star formation occurs at [Fe/H] < —2.5. In 
the scenario of monolithic collapse this occurs 
because of insufficient rate of star formation 
at the epoch of initial contraction due to low 
density of the protogalactic cloud. 

In the computations of Galaxy formation 
according to the monolithic scenario, i.e. in 
model A (Fig. 1), the reason for cessation of 
star formation is the increase of temperature, 
which is due both to gas compression during 
the collapse and to the viscous dissipation. 
Formation of a shock wave as a result of col- 
lapse or of merging of galaxies may serve as a 
physical reason for increase of viscous dissipa- 
tion. Larson [43] has shown that the merging 
of galaxies results in the suppression of gas 
cooling and that this may cause cessation of 
star formation. In the hierarchical scenario 
formation of a galaxy may occur as a result 
of a merging of lower-mass galaxies. Thus, 
an elliptical galaxy with mass > 10^^ Mq 
may have two or three significant (i.e., similar 
mass) precursors and the epoch of its "assem- 
bly" corresponds to the redshift ^ ~ 1 [44], 
i.e., to the stellar age of ~ 9 Gyr. Most of 
the precursors are formed at z ~ 2 — 3 [44] 
experiencing at this time a burst of star for- 
mation, during which the abundances of a- 
elements and heavy elements increase. The 
mergings may be important for disk galaxies 
too, but in this case the mass of the consumed 



galaxy has to be much less, since otherwise 
an elliptical galaxy will form. The metallici- 
ties close to solar one that are discovered for 
some thick disk stars [45] also favor hierar- 
chical scenario with independent evolution of 
precursors. More, it is possible to show that 
the stars of halo and thick disk formed si- 
multaneously [42] and formation of halo stars 
took ~ 1.5 X 10^ yr [46]; the last circumstance 
also does not contradict hierarchical scenario. 

In the hierarchical scenario, the history of 
star formation for our Galaxy that is con- 
sistent with observations may result from 
a merging with a satellite with mass ~ 
2 X 10^° Mq. Figure 2 shows results of 
computations of the evolution of a satellite 
galaxy with the mass of 2 X 10^° Mq, radius 
6.32 kpc, and initial semi-thickness 6.32 kpc 
(these parameters follow TuUy-Fisher rela- 
tion M/R^ = const, where M and it! are 
taken from model A); we designate this model 
as B. In model B beginning of computa- 
tions corresponds to the Universe age 1.5 Gyr 
{z = 2.66). It is seen that after the phase 
of the initial collapse a burst of star forma- 
tion occurs which leads to the increase of oxy- 
gen and iron abundances to [0/H] ^ —0.25 
and [Fe/H] pa —0.5 when the galaxy becomes 
1.5 Gyr old. One may infer that merging 
with satellite occurred at z = 1.5 (the age 
of the Universe 3.5 Gyr). Merging was mod- 
eled within scenario of monolithic collapse 
(model G) using the following initial condi- 
tions: mass of the galaxy 2.2 x 10^^ M©, ra- 
dius - 20 kpc, age of the Universe - 3.5 Gyr. 
At the merging moment, abundance of oxy- 
gen in model A was [0/H] pa —2.2, abun- 
dance of iron - [Fc/H] « —2.6 (Fig. 1), 
while in model B these abundances were 
[0/H] ^ -0.25 and [Fe/H] ^ -0.5, respec- 
tively (Fig. 2). Gas density in our Galaxy 
was 0.25 Mq /pc^, while in the satellite it was 
0.45 Mq /pc^. This means that in the galactic 
region where the merging occurred, mixing of 
the gas would result in abundances of oxygen 
and iron corresponding to the abundances in 
the sateUite galaxy. This was the reason to 
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FIG. 4: Evolution of the chemical composition of the Galaxy in the hierarchical scenario. We show 
together the plots of [0/H] — [Fe/0] from Figs. 2 and 3. The arrow indicates location of the pause 
in star formation. Halo star are shown by squares, thick disk stars by circles, and thin disk stars 
by triangles. Dashed line shows results of computations for the standard model [12]. 



set initial chemical composition of the gas 
similar to the composition in the model B 
at the end of computations. 

Figure 4 shows relation [0/H] — [Fe/0] 
composed from the plots of evolution of abun- 
dances of oxygen and iron in models B and 
C, as well as positions of halo, thick and thin 
disk stars according to Gratton [39]. Loca- 
tion of the pause in star formation is indi- 
cated by the arrow. As Fig. 3 shows, the 
pause in star formation corresponds to the 
range of stellar ages 9 — 10.5 Gyr, in agree- 
ment with observations. 

IV. CONCLUSION 

It is evident that in the star formation 
law it is necessary to account for the turbu- 
lent energy that results from the turbuliza- 
tion of the ISM by supernovae. It is also ev- 
ident that the parameters of the ISM (den- 
sity, turbulent energy, temperature) are sig- 
nificantly different at different length scales. 
In the present study a star formation func- 
tion that depends on the turbulent energy of 



interstellar medium is suggested. The model 
is based on application of Jeans criterion for 
power-law distribution of density perturba- 
tions in the ISM. Such an approach results in 
the Schmidt-type dependence of the SFR on 
the gas density with power ^ 2 and an inverse 
quadratic dependence on turbulent energy of 
the ISM (see Eq. (10)). Applying Jeans cri- 
terion, the model may be generalized to the 
case when a significant infiuence of the chem- 
ical composition of the ISM, magnetic field 
or galactic rotation are assumed. As well, 
we constructed a model for the energy dissi- 
pation for the case when the ISM is consid- 
ered as turbulent over different length scales. 
Star formation and dissipation laws we sug- 
gest, are intended for application in the nu- 
merical models of galactic evolution. Within 
single-zone galactic evolution model we have 
explained the delay of star formation in the 
stellar ages range from 8 — 9 to 10 — 12 Gyr. 
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